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A fully microscopical simulation of nuclear collisions by a new QMD model 
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Abstract. Nucleon-ion and ion-ion collisions at non relativistic bombarding energies can be described by means 
of Monte Carlo approaches, such as those based on the Quantum Molecular Dynamics (QMD) model. We have 
developed a QMD code, to simulate the fast stage of heavy-ion reactions, and we have coupled it to the de-excitation 
module available in the FLUKA Monte Carlo transport and interaction code. The results presented in this work span 
the projectile bombarding energy range within 200 - 600 MeV/A, allowing to investigate the capabilities and limits of 
our non-relativistic QMD approach. 



1 General framework and motivation of this work 
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Fig. 1. 



Double-differential neutron production cross-section 
for Ar projectiles impinging on C at 560 MeV/A bombarding 
energy. The results of the theoretical simulations made by 
QMD + FLUKA de-excitation are shown by solid histograms, 
while the experimental data taken from ref. [14| are shown 
by circles. Distributions at 5°, 10°, 20°, 30°, 40°, 60° and 80° 
(lab) angles are scaled by decreasing even powers of 10. 



expansion stage, characterized by the formation of hot excited 
fragments (pre-fragments). The whole phase occurs on a time 
scale within a few hundreds fm/c, depending on the size of the 
colliding systems and the bombarding energy, and is called the 
"fast" stage of the reaction. 
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Fig. 2. Double-differential proton production cross-section 
for p projectiles impinging on C at 300 MeV/A bombarding 
energy. The results of the theoretical simulations are shown by 
solid histograms, while the experimental data from ref. lfT31l 
are shown by squares. Distributions at 20°, 30°, 40°, 50°, 60°, 
75°, 90° and 105° (lab) angles are scaled by 10 4 , 10 3 , 10 2 , 10 1 , 
10°, 10 -1 , 10~ 2 , 10" 3 , respectively. 




A fully microscopical simulation of nucleon-ion and ion- 
ion collisions, at nucleon level, can be performed, among 
several different approaches |fl~), by means of Quantum Mo- 
lecular Dynamics (QMD) models [2|. They are dynamical 
models which allow to study the phase-space evolution of the 
projectile-target colliding systems, from their initial mutual 
trajectory influence and eventually their overlap, depending on 
the impact parameter, to the compression phase, accompanied 
by a temperature and density increase, up to the following 
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Additionally, improved versions of QMD models (e.g. the 
CoMD one, developed by Papa et al. |3|) allow to compute 
the system evolution even for a longer time, up to thousands 
fm/c, and thus have also been used to describe pre-fragment 
de-excitation, at least in its initial stage. This has led to direct 
comparisons of the results of improved QMD simulations to 
experimental data concerning fragment emission distributions. 
On the other hand, pre-fragment de-excitation can occur on 
a time scale even larger (up to ~ 10~ 15 s). Thus, a complete 
treatment of this slow stage can be covered by different 
models, generally based on statistical considerations. The 
underlying assumption in applying one of these statistical 
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Fig. 3. The same as fig. [2] for p impinging on C at 392 MeV. 
Distributions at 20°, 25°, 40°, 50°, 75°, 90° and 105° (lab) 
angles are scaled by 10 3 , 10 2 , 10 1 , 10°, 10" 1 , 1(T 2 , 1(T 3 , 
respectively. 
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models to nuclear systems is that they are thermalized. While 
at the lowest energies the colliding ions stay close to each 
other for a time long enough for thermalization to occur, to 
define a temperature for the whole system at higher energies 
(several tens MeV/A) can be very problematic, since the 
expansion phase can begin before a global thermalization 
process is completed. Anyway, in the last case, at advanced 
time in the expansion stage pre-fragments are well separated 
and each of them is supposed to be thermalized. Whereas 
theoretical models allow to compute an excitation energy for 
each pre-fragment, from the experimental point of view the 
problem of the determination of hot fragment temperatures 
is still open. On the other hand, planned applications (such 
as hadrontherapy and space radioprotection) need models and 
tools to calculate doses to human bodies and equipment, 
due to radiation exposure. In particular, reliable calculations 
of physical doses require an accurate description of nuclear 
interactions. Nuclear reaction models used for predictions in 
these applications, should be capable of reproducing available 
experimental data concerning particle and fragment emission. 



This paper presents the results of simulations performed 
with a newly developed QMD code, coupled to the de- 
excitation module available in one of the most widely used 
Monte Carlo transport and interaction codes, the FLUKA 
code II4I51 . In particular, FLUKA includes algorithms al- 
lowing to compute evaporation, fission, fragmentation and 
Fermi break-up (for light nuclei) followed by y emission, 
for whichever excited pre-fragment. Those algorithms are 
part of the FLUKA general hadron-nucleus (and lepton- 
nucleus) interaction model called PEANUT (PreEquilibrium 
Approach to NUclear Thermalization) [6|. For further infor- 
mation, the interested reader can check also refs. B7I8I9I101 
and the FLUKA website http://www.fluka.org The excited pre- 
fragments computed by the QMD model at the end of the fast 
stage are then passed through these algorithms for the final 
de-excitation and fragment generation. 
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Fig. 4. Charged projectile-like fragment production cross- 
sections for Si ions impinging on C at 266 and 344 MeV/A 
bombarding energies, respectively in the upper and in the 
lower panel: theoretical predictions are compared to exper- 
imental data from ref. 1161 . Detector acceptance has been 
included in the theoretical simulations. 



2 Comparisons with experimental data 



2.1 Neutron emission 



As far as neutron emission is concerned, performances of 
the model in predicting double-differential neutron production 
cross-sections in the projectile energy range from several tens 
MeV/A up to ~ 400 MeV/A have already been shown in 
refs. Ml 1I12I131 . The QMD we have developed so far includes 
non-relativistic potentials, thus it is expected to work at ener- 
gies up to a few hundreds MeV/A, limited also by the inability 
to deal with pion production and reinteraction. In the effort of 
better understanding up to which energies it can be reasonably 
applied, and which modifications have to be included to 
improve it, a few simulations at energies above ~ 500 MeV/A 
have also been performed. As an example, the neutron double- 
differential production cross section for Ar + C at 560 MeV/A 
bombarding energy is shown in fig. Q] The experimental data 
taken from ref. |[T4l are plotted together with the theoretical 
curves. It is apparent that, as far as forward emission angles 
are concerned (5 - 10°), the model underestimates neutron 
emission tails from ~ 800 MeV/A up to the highest energies. 
On the other hand, at > 20° emission angles, the agreement 
of the theoretical results with the experimental data is quite 
encouraging, and at angles > 40° the neutron tails are nicely 
reproduced. 
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Fig. 5. The same as fig H for Si ions impinging on Cu at 274, 
344 and 545 MeV/A bombarding energies, respectively in the 
upper, intermediate and lower panel. 
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Fig. 6. The same as figs|4]and[5]for Si ions impinging on Pb at 
274, 364 and 540 MeV/A bombarding energies, respectively 
in the upper, intermediate and lower panel. 



2.2 Proton emission 

While systematic data concerning neutron emission in heavy- 
ion collisions exist, as far as proton emission in ion-ion colli- 
sions is concerned, less data are available in literature. We thus 
start considering nucleon induced reactions, for which many 
more data concerning proton emission exist. Indication of the 
performances of our code in reproducing the experimental 
data presented by ifTBI is shown in fig. [2] and [3] for p + C 
at 300 and 392 MeV bombarding energies, respectively. The 
linear scale on the energy axis allows to better appreciate the 
position of the broad emission peaks. With the exception of 
the quasi-elastic peak, the general agreement is quite good 
as far as forward emission angles are concerned, while the 
neutron tails at larger angles are underestimated, especially 
at the highest energies. The slight departure from smoothness 
of some of the lines obtained by simulation is not due to 
statistics. It is instead an artifact introduced by the fact that 
quite a few QMD ion initial states have been used to perform 



these simulations, where projectiles are single protons. In fact, 
only seven 12 C initial states have been used to obtain these 
figures. As a general result, when nucleon - ion or light ion - 
light ion simulations are performed, better results are obtained 
considering a larger number of QMD initial states. 



2.3 Projectile-like fragment emission 

As far as fragment emission is concerned, a few results 
obtained by our QMD coupled to FLUKA de-excitation have 
already been presented in refs. H 1 21 1 31 . 

In the present paper, the focus is given to the case of 
projectile-like fragments produced by Si ions impinging on 
targets made of C, Cu and Pb. Systematic studies of Si 
fragmentation on targets of light, intermediate and heavy 
composition have been carried out in the last few years at 
the HIMAC in Chiba, and at the BNL AGS. The results of 
these experiments were recently published ifToll . Comparisons 
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Fig. 7. Scatterplot of fragment production cross-sections pre- 
dicted by QMD + FLUKA de-excitation (vertical axis) vs. 
measured cross-sections (horizontal axis) for Si ions imping- 
ing on Pb at 274 (crosses), 364 (squares) and 540 (circles) 
MeV/A bombarding enegies. The errorbars on the theoret- 
ical results are statistical only. The two lines respectively 
above and below the diagonal correspond to deviations of 
the modelled results from the experimental ones by ± 33 %. 
The largest discrepancies are seen for the Z = 13 fragment 
production cross-sections. 



with the predictions of our model, as far as projectile-like 
fragment production cross-sections are concerned, are under 
way. A few results, for the lowest energies cases, are plotted 
in figs. |U |5] IS] Detector acceptance was taken into account 
in performing our simulations by means of cuts in energy 
and angle, which allow to select only projectile-like fragments 
escaping the target within a few degrees around forward 
direction, according to the angular acceptances given in ifTBI . 
Due to the features of the detectors used in these experiments, 
data are available for fragments with charge Z = 6 - 13 only. 
For lower charges (Z < 5), it was very difficult to disentangle 
the peaks corresponding to the separate contribution of each 
fragment species to the energy released in the Si detectors 
used. Sometime, also the broad peaks corresponding to the 
Z = 6, 7 charges were determined with some uncertainties, 
especially at the lowest energies. Anyway, the authors of 
ref. 06 1 observe that their largest uncertainties concern the 
determination of the Z = 13 fragment production cross- 
sections, due to superposition effects with the high-energy 
tails of primary ions that cross the targets without interacting, 
which have to be cut. These uncertainties are more important 
at lower energies, due to the use of thin targets to avoid 
significant energy losses of projectiles, which would lead to 
misleading reaction cross-sections. 

In all cases discussed in the present paper, the theoreti- 
cal model systematically underestimates the Z = 13 cross- 
sections, while for Cu and Pb targets it overestimates the 
Z = 6 production cross-sections. The best agreement with 
the experimental data is instead observed for fragments with 
Z = 10, 11, and for fragments with Z = 12 at energies below 
400 MeV/A. At higher energies, an underestimation of the 
Z = 12 fragment abundances occurs as well. This could be an 
indication that the QMD developed so far gives more reliable 
results when used in the study of central collisions, while it 
is more difficult to apply it to the study of stripping reactions, 
where only one or two nucleons of the projectile ion interact 



with the target. Furthermore, the authors of ref. IfTBI suggest 
the hyphotesis that electromagnetic dissociation can contribute 
to the Z = 12, 13 fragment production cross-sections at high 
energy for the heaviest targets. At present, we have not tested 
this hypothesis yet. Finally, as far as Fluorine fragments 
(Z = 9) are concerned, in all cases the theoretical model 
confirms their suppression, at least from a qualitative point of 
view, with respect to the abundance of the other species close 
in charge, as observed in the experiment. 

To obtain a global overview of the deviations of the 
predictions of our theoretical model from the experiments, 
the results can be summarized in scatterplots, such as the 
one shown in fig [7] for the Si + Pb case, including data 
at different energies for fragments with Z = 6 - 13. With 
the exception of the cases of the heaviest (and lightest) 
charged fragments already discussed, the figure does not show 
particular systematics. For the heavy target considered, it is 
apparent that the discrepancies between the model predictions 
and the experimental data do not exceed ~ 30 - 35% in most 
cases. Theoretical fragmentation models based on QMD codes 
provide, in general, more reliable predictions in case of inter- 
mediate and heavy mass targets than in case of lighter ones. 

Collaboration with F. Ballarini, G. Battistoni, F. Cerutti, A. Fasso, 
A. Ferrari, E. Gadioli, A. Ottolenghi, L.S. Pinsky, J. Ranft and 
PR. Sala is acknowledeged. This work was supported by the Uni- 
versity of Milano. 
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